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Abstract 

Recent experiments have reported measurements of rainbow scattering features in the angular 
distributions of hyperthermal Ar colliding with LiF(OOl) [Kondo et ai, J. Chem. Phys. 122, 
244713 (2005)]. A semiclassical theory of atom-surface collisions recently developed by the authors 
that includes multiphonon energy transfers is used to explain the temperature dependence of the 
measured scattered angular distribution spectra. 
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I. INTRODUCTION 


Recently, Kondo et al. published an extensive series of experimental results on the scatter¬ 
ing of well-defined incident beams of liyperthermal Ar atoms from clean LiF(OOl) surfaces.- 1 ^ 
The data consist of angular distributions as a function of final scattering angle for an exper¬ 
imental apparatus having a fixed angle of 90° between the incident beam and the detector 
position, and at each final angle all scattered particles were detected regardless of trans¬ 
lational energy. In contrast to many other measurements involving liyperthermal surface 
scattering of heavy rare gas atoms,- these measurements are noteworthy because they were 
carried out for a large range of incident energies, from 315 to 705 meV, and temperatures 
ranging from 300 to 573 K. Under many conditions, and particularly for scattering in the 
plane of the (100) direction, the measured angular distributions exhibited two broad, but 
quite distinct peaks attributed to classical rainbow scattering caused by the corrugation of 
the LiF(001) surface. 

In this paper we wish to analyze the Ar/LiF(001) data using using a new semiclassical 
atom-surface scattering theory that has been demonstrated to explain experimental mea¬ 
surements for liyperthermal energy scattering of rare gas atoms from liquid metal surfaces.-^ 
This theory, briefly described in somewhat more detail in Sec. [TT] below, uses an interaction 
potential that is relatively straightforward and simple, a corrugated hard repulsive wall that 
allows for phonon vibrations due to the thermal motions of the underlying surface atoms. 
Thus, although the interaction potential as a function of the ^-coordinate normal to the sur¬ 
face is simple, the surface periodicity is retained and the the exchange of multiple phonons 
in the classical scattering limit is treated realistically. The results of our calculations explain 
the measured temperature dependence of the angular distributions, and show that at fixed 
incident energy, but with increasing surface temperature the increasing inelastic background 
broadens the rainbow peaks and makes them less distinct. 

As a function of increasing incident energy it was noticed that the angular separation 
between the two rainbow peaks increased, giving the appearance of a surface corrugation 
height that decreased with energy.- However, it was predicted some years ago that such be¬ 
havior of rainbow features could be the result of refraction of the incident and final scattering 
beams by the attractive adsorption potential well in front of the repulsive barrier.- 

The interesting energy dependence of these data has been analyzed by Poliak et al. 
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using a classical atom surface scattering perturbation theory that includes energy transfer 
to the phonons via a Langevin thermal bath and with the surface periodicity described 
by a one-dimensional corrugation function that contains multiple Fourier components.- 1 ^ 
By including the Langevin description of energy transfer between projectile and surface, 
they showed that the Ar/LiF(001) angular distributions could be quite well replicated by 
a single, fixed, one-dimensional surface corrugation potential that included an attractive 
adsorption well and a very strong second order sinusoidal Fourier component. Thus, they 
showed that the decrease with increasing incident energy in the angular separation of the 
two rainbow peaks could be explained by a single, fixed corrugation function, and did not 
imply a corrugation that changed with incident energy, ffowever, their corrugation function 
contained not only a sinusoidal term describing one-dimensional periodicity with the LiF 
lattice constant, but also an even stronger second harmonic term with periodicity half that 
of the lattice constant. The use of a corrugation function with strong higher-order Fourier 
components was originally proposed by Tully, with the explanation that the higher order 
terms mimicked the dynamical behavior of the interaction potential, i.e., they helped to 
account for the fact that the interaction potential can vary significantly depending on the 
height above the surface where the incoming projectile is deflected from the surface within 
the periodic unit cell. 9 - In a subsequent extension of this theory to a second order classical 
scattering perturbation treatment Poliak et al. have demonstrated that observed decrease 
in rainbow peak separation with increasing incident energy can be explained with a one 
dimensional corrugation containing only a single Fourier component, but is dependent on 
the interaction potential having an appropriate attractive adsorption well. 10 A very general 
molecular dynamics approach which includes fully three-dimensional scattering has recently 
been reported.— Although at its present state of implementation calculations with this 
treatment shows only qualitative agreement with the Kondo data of Refs, fi] and p], they 
clearly demonstrate the shifts of the rainbow peaks with increasing incident energy. 

In this paper we use a potential that is purely repulsive, without an attractive well, and 
although our calculations can match the energy dependence of the Ar/LiF(001) data, it 
requires a corrugation that decreases in height as the incident energy is increased. In the 
light of Refs. discussed above, this is clearly caused by the absence of an attractive 

well in the potential. However, the temperature dependence of these data, which was not 
examined in Refs, is well explained by our theory of energy exchange with a repulsive 
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barrier. Also presented are several calculations for energy-resolved scattering spectra that 
illustrate how more detailed information on the interaction potential can be obtained from 
such measurements, as opposed to that provided by measurements of the more restrictive 
angular distributions. 


II. THEORY 


The theory used in this paper is fully developed in Ref. f£], so we present only a few salient 
details here. The starting point is the generalized Fermi golden rule for the transition rate 
of a projectile with initial momentum p; going to momentum state p f after the collision 

MP7, Pi) = FT ( \ T M 2 ^ S f ”£) / ’ W 
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where Tfi is the matrix element of the transition operator with respect to final and initial 
states of the system of projectile plus target, E % and Ef are the initial and final energies of 
the entire system, h is Planck’s constant, the angular brackets signify an average over all 
initial states of the surface and the Y2{ nf } indicates a sum over all final states of the target. 

The experimentally measured quantity is the differential reflection coefficient which gives 
the fraction of scattered intensity deflected into the small final energy interval dEf and small 
final solid angle dflj and is related to the transition rate by 
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where on the right hand side the factor pj comes from the Jacobian relating phase space to 
energy space, and the perpendicular component of the incident momentum p iz is proportional 
to the incident particle flux. 

Applying the generalized Fermi golden rule of Eq. (JTJ) to scattering from a vibrating hard 
corrugated wall within the eikonal approximation of Levi et al .—, under conditions where 
large numbers of phonons are excited, leads to a transition rate expressed as a double integral 
over the surface of the target 
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In Eq. 05 the three-dimensional recoil energy AEq is related to the momentum vector 
Pf — pj through AE 0 = (p f — p. ; ) 2 /2Mc where Me is the effective mass of the target, P 
is the component of the momentum vector parallel to the surface, and £(R) is the surface 
corrugation function. The mass of the projectile is m, T$ is the surface temperature, k B is 
Boltzmann’s constant, and v B is a velocity parameter—arising from the law of conservation 
of parallel momentum for all phonons transferred. The proportionality constant C involves 
only factors of the quantum confinement length. 

The effective mass Me is a measure of how many LiF molecules are involved in the recoil 
of the surface as a result of the collision with the Ar projectile. The parameter vr is a 
weighted average over phonon speeds parallel to the surface^h^ and its significance, as well 
as suggestions of how to best measure it, are discussed in Refs. pi - — ]. 

Eq. (J3]) is a semiclassical approximation, but evaluated in the limit in which large numbers 
of phonons are transferred in the collision, so many phonons that the specific details of 
the phonon distribution function become unimportant and the energy transfer is expressed 
through the recoil function AEq. For large projectile energies the coherence length for 
quantum interference becomes negligibly small, and the double integral over the surface can 
be replaced by integrals over a single unit cell, multiplied by the number of unit cells which 
is simply a large constant. 

An interesting limiting case arises upon assuming that the surface is flat, which means 
the corrugation function £(R) is constant. Taking the classical limit in this case the integrals 
can be carried out leading to 
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With the exception of the factor of p 2 z Eq. 05) is the form first derived by Brako and 
Newns . 13 > 14,18 ~ — The original Brako-Newns theory neglected the effects of the elastic, or 
non-vibrating, part of the interaction potential. Later work has shown that the elastic part 
of the potential introduces a multiplicative factor which is a squared matrix element of the 
elastic potential. The factor of pj z comes from the this matrix element evaluated for a hard 
repulsive wall as introduced by the eikonal approximation. Eq. 05) has been successfully 
used to describe the energy and temperature dependence of atom-surface scattering for 
many systems observed under classical conditions.— —^ However, because it describes the 
surface as a flat repulsive wall with no static corrugation, it is incapable of explaining rainbow 
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features. The theory embodied in Eq. (]3j) is a generalization of the Brako-Newns expression 
gD to include surface corrugation. 


III. COMPARISON WITH EXPERIMENTS 

For the calculations reported here we use a one dimensional model for the surface corru¬ 
gation with a sinusoidal corrugation function 

£(R) —> f(z) = hbcos(^^j . (5) 

The lattice parameter of LiF is a = 4.02 A->22 thus calculations based on Eq. (J3J) depend 
on three parameters. These are the dimensionless corrugation parameter h, the effective 
mass Me of the LiF, and the velocity parameter Vr. It is also noted that even though the 
corrugation is one-dimensional the scattering theory of Eq. (|3]) is fully three-dimensional 
and energy exchange involves all three directions on the surface. 

Angular distribution spectra for in-plane scattering of Ar from LiF(OOl) along the (100) 
direction are compared with calculations as a function of final detector angle Of in Fig. |U 
The spectra were measured at three different surface temperatures of 300, 435 and 573 K 
as marked and the incident energy is E % = 525 meV. At the lowest temperature of 300 K 
there are two clearly distinct rainbow peaks, asymmetrically positioned about the specular 
position Of = 7t/4, and those two peaks differ significantly in intensity. As the temperature 
is raised the increased inelastic scattering causes the two-peak feature to become less distinct 
and the whole scattered distribution broadens slightly. 

It is of interest to briefly discuss the expected rainbow pattern from an elastic, mirror- 
reflecting surface. For scattering in the plane of incidence with a symmetrically shaped 
corrugation, such as the cosine function of Eq. (J5]), the rainbow positions would be sym¬ 
metrically positioned about the specular position. If the angle of incidence 0, is fixed, the 
rainbow angles corresponding to mirror reflection from the inflection regions of the cosine 
corrugation of Eq. (J5]) would appear at the two final angles Of = 0* ± 2 tan^ 1 (27r/r). For the 
the fixed source-detector geometry in the experiments considered here the mirror reflection 
rainbow angles of the cosine corrugation would appear symmetrically about the specular 
position with Of = 7t/4 ± tan^ 1 (2'/rh). Clearly the data shown in Fig. |T]do not exhibit the 
symmetry of a mirror surface. The subspecular rainbow peak appears much closer to the 
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specular position Of = 7r/4 than the supraspecular peak. The large inelastic energy exchange 
between the Ar and the LiF surface causes the rainbow pattern to appear quite differently 
from the form predicted by mirror reflection. The softness present in an interaction potential 
more realistic than the present hard wall repulsion has also been demonstrated to contribute 
to the shifts of the rainbow peaks,7’ 8 ’ 10,11 but in the present calculation all of the shift is 
caused by inelastic energy transfer. 

The calculations from Eq. d3J) are shown as solid curves in Fig. [0 These are obtained 
from the differential reflection coefficient of Eq. (J2]) after integrating over all final energies. 
In the integral over final energies a factor of 1 /\[E] is introduced to account for the 1 jvf 
velocity dependence of the experiment which uses a density detector. The calculations use 
the following parameter values: h = 0.024, Me is 13 times the 25.9 arnu mass of LiF or 
337 amu, and v r = 2000 m/s. These same mass and vr values are used for all calculations 
shown in this paper. The calculated intensities are not strongly affected by the choices of 
these latter two parameters. If Me is changed by 20% and/or vr is changed by 50% the 
results differ very little, hence the results are not critically dependent on the fitting of these 
two parameters. The repeat distance for the one-dimensional corrugation was chosen as 
b = a/2 = 2.01 A. This choice is motivated by the fact that for the LiF(001) surface the 
predominant feature causing the corrugation is the bumps presented by the large F ions, 
which form a square array, with sides of length a/\/ 2 = 2.84 A, rotated 45° with respect 
to the LiF(001) unit cell. He atom diffraction peaks in the (100) direction are separated by 
parallel wave vectors of 3.1 A -1 which corresponds to diffraction from rows with a spacing 
of b = a/2. 1 ' 22 ' 23 Nevertheless, this choice is not important in the calculations presented here 
for the following reason: as mentioned above for the case of classical scattering from a hard 
mirror surface the positions of the rainbow peaks are defined by the slope of the corrugation 
at its inflection point, and from Eq. (J5]) this is d^/dx = ±2irh and is independent of the 
repeat distance b, but depends on the dimensionless corrugation parameter h. Because 
the present calculations are in the classical limit where quantum interference is negligible 
our calculations are similarly sensitive primarily to the corrugation parameter h. Thus 
calculations using a repeat distance b = a are nearly identical to those for b = a/2. 

The calculations are strongly sensitive to the choice of the dimensionless corrugation 
parameter, as will become evident in the discussion of below. The value of the corrugation 
parameter h = 0.024 is chosen to give the best fit with the data at the lowest temperature 
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of 300 K. Because the data were reported in arbitrary units, the calculations are normalized 
to the data at a single point in each panel of Fig. [0 At all temperatures the calculations 
match the data well and appear to explain the smoothing and broadening of the rainbow 
features as the temperature increases. The asymmetry in both angular position and relative 
intensity of the two rainbow peaks is well predicted, and the temperature broadening is also 
reasonably explained by the inelastic transfers built into this theory. 

Fig. [2] shows the temperature dependence of the angular distribution spectra measured 
in the (110) direction for E t = 525 meV. The (110) direction lies along the squares formed 
by the fluorine atoms on LiF(001) so the one dimensional sinusoidal period is chosen to be 
b = a/y /2 = 2.84 A. The three surface temperatures are the same as in Fig. [D In this 
direction the expected two peaks of the rainbow pattern have merged together and appear 
as a single broad feature, an indication that the apparent corrugation strength is weaker. 
The calculations are similar to those of Fig. |Tj except now the corrugation parameter takes 
the smaller value h = 0.01. With this corrugation the temperature dependence of the data is 
explained rather well which indicates that the merging of the two expected rainbow peaks and 
the observed broadening of the single merged feature may be explained as a consequence of a 
rather large multiphonon energy transfer from the Ar to the surface. This statement should 
be tempered with caution, however, because it has been shown that softness in the repulsion 
for a more realistic interaction potential, as well as the contribution of the attractive potential 
well, can also contribute to the shifts and broadening of these features.--— The calculations 
also indicate that the apparent corrugation in the (110) direction is significantly smaller 
than in the (100) direction. 

Although it is of interest to examine the temperature dependence of angular distribu¬ 
tions, the following figures present calculations that show how substantially more detailed 
information about the collision process can be obtained from energy-resolved measurements. 
Fig. [3] shows the calculated temperature dependence of the scattered intensity as a function 
of final energy for the same conditions as for the angular distributions in Fig. |Tj except that 
the detector position is held at the specular position, i.e., 9,, — 9f — 45°. At all detector 
angles 9f for which the scattered intensity is appreciable the calculated energy-resolved plots 
consist of a single, broad peak. The position of the most probable final energy (the peak 
position) remains the same at all temperatures, but the most probable peak intensity (the 
peak height) decreases strongly with increasing temperature while the full width at half max 
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(FWHM) increases. This is the expected behavior since the calculations are unitary, imply¬ 
ing that the scattered flux must equal the incoming flux. Thus, as the temperature increases 
and distributes the scattered particles over a larger range of final energies the maximum scat¬ 
tered intensity must decrease in order to preserve the total flux. It is also noticed that at 
higher temperatures the scattered intensity develops an increasingly asymmetric tail on the 
high energy side. This indicates that at higher temperatures increasingly larger numbers of 
atoms are “heated up” and some are predicted to leave the surface with substantially larger 
energies than the incident beam, although the stable peak position indicates that the most 
probable energy loss under these conditions remains at about 10% of the incident energy. 
The FWHM increases with temperature at a rate much faster than a proportionality to y/Ts 
that would be expected from examination of the smooth surface model of Eq. (J4j) and this 
is apparently due to the rapidly increasing asymmetry in the high energy tail 

It is of interest to note the sharp decrease in most probable intensity as a function 
of temperature in the energy-resolved spectra shown in Fig. [3l The calculated angular 
distributions, such as shown in Figs. H] and [2j also show a decrease in maximum peak 
intensities, but not nearly so strong as for typical energy-resolved spectra. This is because an 
angular distribution is a measurement that, at a given detector position, sums the intensities 
over all final energies, and again because of the unitarity condition this means that the 
temperature dependence will be less pronounced. 

It has been demonstrated that the temperature dependence of an energy-resolved mea¬ 
surement can be directly related to the surface corrugation amplitude,-^ 24 and this is il¬ 
lustrated in Fig. |4] which shows as open circles the most probable intensity of Fig. [3] as a 
function of Tg. Also shown in Fig. [4] as a solid curve is the expected behavior of approxi¬ 
mately 1/T^ 2 for the smooth surface model as seen from the prefactor to the exponential of 
Eq. ((4]) , and the dashed curve shows the expected behavior of a highly corrugated surface of 

i /n 

discrete scattering centers which goes approximately as 1 /T s . Calculations for a surface 
with corrugation intermediate between these two extremes should give a locus of points on a 
curve lying between the 1/T^ 2 and 1/T X J 2 behaviors, and this is demonstrated by the open 
circle points taken from Fig. [3] The relative position of the curve traced by those points is 
defined by the value of h for the corrugation amplitude of the corrugation in Eq. (J5]) above. 
The proximity of the present calculations to the 1/T^ 2 behavior indicates that this system, 
as expected, is rather highly corrugated. 
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As another example demonstrating the potential sensitivity of energy-resolved spectra 
to the parameters of the interaction potential, Fig. (J5]) shows calculations under the same 
conditions as Fig. [3] at the temperature of 573 K, but with three different effective masses. 
As stated above, the angular distributions of Figs. Q]and[2]are rather insensitive to the value 
chosen for the effective mass, and the calculated angular distributions under those conditions 
would be essentially unchanged if the mass were increased or decreased by as much as 20%. 
However, the situation is quite different for the energy resolved spectra, as seen in Fig. (J5]l 
which shows, in addition to the spectrum for the Me value of 13 LiF masses, curves for 
Me = 11 and 15. Increasing or decreasing Me by ~20% gives rise to most probable final 
energies that differ by as much as 50 meV which is a very readily measurable amount. Thus 
it is very apparent that energy-resolved spectra are a much more sensitive way to sample 
the effective mass of the target. 


IV. DISCUSSION AND CONCLUSIONS 

In this paper we examine a set of data for angular distributions of hyperthermal energy 
Ar atoms scattering from the LiF(OOl) surface. These data are analyzed with a new theory 
of atom scattering from corrugated surfaces under conditions for which the collision is rapid 
and results in the transfer of large amounts of energy via multiphonon excitations. Previ¬ 
ously, we showed that relatively straightforward temperature-dependent measurements on 
the maximum intensities of energy-resolved scattering spectra, taken under conditions of 
fixed incident and detector angles, could be used to estimate the corrugation height of the 
surface.— ^ The presently analyzed measurements are angular distributions, meaning that at 
each angle all particles are detected regardless of final energy. This provides less detailed 
information about the scattering process because of the integration over all final energies, 
but the detection process is simpler and can be followed over a range of final angles and 
thus is capable of measuring effects such as the classical rainbow peaks discussed here. 

The angular distribution data exhibit a number of characteristic effects. Many of the 
angular distributions show two distinct rainbow peaks, the signature of a simply corrugated 
surface with two inflection regions. However, these two peaks are not symmetrically posi¬ 
tioned in angular separation with respect to the specular direction as would be expected 
for reflection from a mirror surface. As the surface temperature is increased the rainbow 
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pattern becomes less distinct and at the highest temperatures, and especially for the (110) 
direction, they merge into a single broad structure. 

Within the confines of our simple interaction potential, the calculations indicate that 
the positional asymmetry of the two rainbow peaks, in both the high-symmetry (100) and 
(110) directions, can be largely accounted for due the inelastic transfers resulting from the 
exchange of energy. This is consistent with previous work of others,- 7,10 but there they found 
that some of the asymmetry is due to softness and the attractive potential well that would 
appear in a more realistic interaction potential. Because our interaction potential lacks 
softness and has no well, it is possible that this could result in over- or underestimation of 
our chosen interaction parameters. 

A significant contribution of the present work is that it explains the temperature depen¬ 
dence of the measured angular distributions. In both symmetry directions, our calculations 
match the spectra reasonably well at low temperatures. As the temperature is increased our 
calculations show how the rainbow peaks become less distinct and gradually merge into a 
single, broad feature as a consequence of larger energy transfer with the surface. 

It should be noted that this work is based on using a one-dimensional corrugation, al¬ 
though the scattering theory and calculations are fully three-dimensional. The real corruga¬ 
tion of the LiF(OOl) surface is two-dimensional, which could have distinct and qualitatively 
different scattering effects as has been indicated by recent work.-! 1 Although this work, as 
well as previous analysis, 7,10 arrives at physically reasonable explanations of the Ar/LiF(001) 
angular distributions, these results need to be verified with a more realistic three-dimensional 
treatment. Work is presently in progress to carry out such three-dimensional calculations 
within the framework of the present theory. 

In addition to direct comparison with the experimentally measured angular distribu¬ 
tions, we showed in Figs. [3]|5] several calculations of predicted behavior of energy-resolved 
spectra taken at fixed incident and detector angles. These examples demonstrate that 
the energy-resolved spectra are substantially more sensitive to experimental parameters 
than are the energy-integrated angular distributions. This indicates that the best types of 
experiments to give the most information on the parameters characterizing the interaction 
are measurements of energy-resolved spectra at each combination of final and initial angles. 
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Temperature dependent angular distribution spectra for Ar scattering from 
LiF(001)(100) with an incident energy of 525 meV and three different sur¬ 
face temperatures as marked. Solid curves are the present calculations with 
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Temperature dependent angular distribution spectra for Ar scattering from 
LiF(001)(110) with an incident energy of 525 meV and h = 0.01. All other 

conditions are the same as in Fig. [I] . 

Calculations showing temperature dependence of energy-resolved spectra for 
Ar scattering from LiF(001)(100) at an incident energy of 525 meV. The 
incident and final angles are fixed at 0j — Of — 45°, while all other conditions 

are the same as for Fig. [0 . 

Temperature dependence of the most probable intensity of the energy-resolved 
spectra shown in Fig. [31 compared to that for smooth and discrete surfaces, 

shown as solid and dashed curves, respectively. 

Calculations of energy-resolved spectra for Ar scattering from LiF(001)(100) 
at T s = 573 K showing dependence on effective surface mass Mq■ Solid curve: 
same as Fig. [3j dashed curve: with effective mass reduced by the mass of 2 
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FIG. 1: Temperature dependent angular distribution spectra for Ar scattering from LiF(001)(100) 
with an incident energy of 525 meV and three different surface temperatures as marked. Solid 
curves are the present calculations with a corrugation parameter h = 0.024 and the experimental 
points are from Ref. £*]. 
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FIG. 2: Temperature dependent angular distribution spectra for Ar scattering from LiF(001)(110) 
with an incident energy of 525 meV and h = 0.01. All other conditions are the same as in Fig. [T| 
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FIG. 3: Calculations showing temperature dependence of energy-resolved spectra for Ar scattering 
from LiF(001)(100) at an incident energy of 525 meV. The incident and final angles are fixed at 
6i = Of = 45°, while all other conditions are the same as for Fig. [1] 
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FIG. 4: Temperature dependence of the most probable intensity of the energy-resolved spectra 
shown in Fig. [3l compared to that for smooth and discrete surfaces, shown as solid and dashed 
curves, respectively. 
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FIG. 5: Calculations of energy-resolved spectra for Ar scattering from LiF(001)(100) at T$ = 573 
K showing dependence on effective surface mass Me- Solid curve: same as Fig. [3l dashed curve: 
with effective mass reduced by the mass of 2 LiF molecules; dotted curve: with mass increased by 
two LiF molecules. 
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